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Abstract 

We show that the Standard Model Lagrangian, including small neutrino masses, 
, has an anomaly-free discrete Zq symmetry. This symmetry can emerge naturally from 

(ij^ + Li + Lj — 2Lfc) gauge symmetry (Lj is the ith lepton number) and ensure the 
stability of the nucleon even when the threshold of new physics A is low. All AS = 1 
and AB = 2 (B is the baryon number) effective operators are forbidden by the Z§ 
symmetry. AB = 3 operators are allowed, but they arise only at dimension 15. We 
estimate the lifetime for "triple nucleon decay" resulting from these operators and 
find that A can be as low as 10 2 GeV. We suggest a simple mechanism for realizing 
Q_i' reasonable neutrino masses and mixings even with such a low scale for A. 
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1 Introduction 



The Standard Model (SM) has been highly successful in explaining all experimental obser- 
vations in the energy regime up to a few hundred GeV. However, it is believed to be an 
effective field theory valid only up to a cutoff scale A. Non-renormalizable operators which 
are gauge invariant but suppressed by appropriate inverse powers of A should then be con- 
sidered in the low energy effective theory lj. The dimension 5 operator ££HH/Al (£ is the 
lepton doublet) which violates lepton number (L) by two units is the lowest dimensional of 
such operators. Experimental evidence for neutrino masses suggests the effective scale of 
L- violation is around A L ~ fO 14 — fO 15 GeV. The d = 6 operator QQQi/ ' K 2 B violates both 
baryon number (B) and lepton number and leads to the decay of the nucleon. The current 
limits on proton lifetime are r v > 5 x fO 33 yrs for p — > e + 7i° and r p > 1.6 x fO 33 yrs for 
p — > vK + |2j. These limits imply that As > fO 15 GeV. Grand Unified Theories with or 
without supersymmetry generate such B- violating operator with A# ~ fO 14 — fO 16 GeV. 
These theories are currently being tested through nucleon decay. Any new physics with a 
threshold A less than the GUT scale will thus be constrained by both proton lifetime and 
neutrino masses [TJ E3 E] ■ 

It is widely anticipated that new physics will show up around the TeV scale. This 
expectation is based on the stability of the Higgs mass. Supersymmetric theories with 
threshold of a few hundred GeV have been extensively studied as a solution to the Higgs 
mass problem (or the hierarchy problem). Indeed, these theories allow for both B and 
L violating operators. Care must be (and usually is) taken to satisfy the experimental 
constraints. 

A second scenario that has been widely discussed recently to address the hierarchy prob- 
lem is the possible existence of large extra dimensions Here the fundamental scale of 
gravity Mpi is lowered from its 4 dimensional value of 10 19 GeV to M* ~ a few TeV. Non- 
renormalizable operators associated with quantum gravity which are only suppressed by 
inverse power of M* will then be present in the low energy theory. L and B violation would 
put major constraints on these theories. Similar remarks will apply to other solutions to the 
hierarchy problem, such as composite Higgs, little Higgs 6J, etc. 

In this paper we address the role discrete gauge symmetries [7] can play in suppressing 
baryon number and lepton number violating processes even when the scale of new physics is 
low. The SM effective lagrangian does not have a continuous anomaly-free symmetry that 
can suppress these processes. This is our reason for focusing on discrete symmetries. It is 
preferable that such symmetries have a gauge origin since all global symmetries are expected 
to be violated by the quantum gravitational effects [S]. We find that the SM lagrangian, 
including small neutrino masses, has an anomaly-free Z§ symmetry. This Z§ acts as discrete 
baryon number and forbids potentially dangerous AB = 1 and AB = 2 nucleon decay 
processes. AB = 3 operators are allowed, but they arise only as d — 15 operators. "Triple 
nucleon decay" can proceed through these operators; we estimate its rate and find A B can 
be as low as 10 2 GeV. We show that the Z 6 symmetry has a natural gauge origin in the 
anomaly-free U(l) corresponding to Ir + L{ + Lj — where Li is ith lepton number and 
In is the third component of the righthanded isospin. No new particles are introduced to 
cancel anomalies except the righthanded neutrinos. We show how an additional Z5 or Zj 
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gauge symmetry generates the desired neutrino masses and mixings even with a low scale 
for A. 

Discrete gauge symmetries have been utilized in suppressing nucleon decay as well 
as in addressing other aspects of physics such as solving the fi problem of supersymmetry, 
fermion mass hierarchy problem [TTH ITT] and the stability of the axion ^JJ. A Z 3 baryon 
parity was found in Ref. that suppresses nucleon decay. In order for it to have a gauge 
origin, extra particles were introduced. 4 There have been also attempts to gauge baryon 
number by introducing a fourth family of fermions or with the use of extra dimensions 
[T3] . The spirit of this paper is similar to that of Ref. [UJ C21 EE3j , our models are however 
much more economical. 



2 Gauged Baryon Parity and Suppression of Proton 
Decay 

In this section we show that the SM lagrangian has a discrete Z e gauge symmetry which 
forbids all AB = 1 and AB = 2 baryon violating effective operators. This can be seen as 
follows. The SM Yukawa couplings incorporating the seesaw mechanism to generate small 
neutrino masses is 

C Y = Qu c H + Qd c H* +£e c H* + £u c H + M R u c u c . (1) 

Here we have used the standard (lefthanded) notation for the fermion fields and have not 
displayed the Yukawa couplings or the generation indices. This lagrangian respects a Zq 
discrete symmetry with the charge assignment as shown in Table ^ Also shown in Table 
^are the charge assignments under the Z3 and Z2 subgroups of Zq. The Z% assignment is 
identical to that in Ref. jU] 
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Table 1: Family-independent Zq charge assignment of the SM fermions and the Higgs boson along with 
the charges under Z3 and Zi subgroups. 



From Table Q it is easy to calculate the Z§ crossed anomaly coefficients with the SM 
gauge groups. We find the SU(3)c or SU(2)l anomalies to be 

A [SU(3)c] 2 xZ 6 = 3N g 

A{su(2) L } 2 xz 6 — Ng (2) 

4 The charges of the SM particles under the U(l) symmetry in Ref. [5] are {Q, u , d, I, e} = 
(0, — 1, 1, — 1, — 1). The model has extra leptons T(3), l'(0), e'(3), e'(0) (where the U(l) charges are given 
in parentheses) as well as exotic leptons with electric charge of ±1/3 and 0: 2 x F 1 / 3 (3), 1 x F _1 / 3 (3), 
2 x F-V 3 (0), 1 x F 1 / 3 (0) and 5 x iV°(-3), iV°(l), N°(3). 
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where N g is the number of generations. The condition for a Z^ discrete group to be anomaly- 
free is 

N 

Ai = — mod N (3) 

where i stands for SU(3)c and SU(2)l- For Z 6 , this condition reduces to Ai = 3 mod 6, so 
when N g = 3, Z e is anomaly-free. Obviously, the Z 3 and Z 2 subgroups are also anomaly- 
free. The significance of this result is that unknown quantum gravitational effects will respect 
this Zq. It is this feature that we utilize to stabilize the nucleon. Absence of anomalies also 
suggests that the Z 6 may have a simple gauge origin. 

To see how the Z$ forbids AB = 1 and AB = 2 processes, we note that it is a subgroup 
of U(1)2y-b+3L where Y is SM hypercharge jUIj. We list in Table El the charges under the 
three U(l) symmetries. Comparing the last line of Table El with the charge assignment of 
Table [T] it is clear that the Z 6 is a subgroup of U(l)zy-B+3L- Any Z 6 invariant effective 
operator must then satisfy 

2 AY - AB + 3AL = mod 6. (4) 

Invariance under U(l)y implies AY = 0. Consider AB = 1 effective operators which must 
then obey (from Eq. (@J)) 3AL = 1 mod 6. This has no solution, since 3AL = mod 3 
from Table El Similarly, AB = 2 operators must obey 3AL = 2 mod 6 which also has no 
solution. AB = 3 operators, which corresponds to 3AL = mod 6, are allowed by this Z 6 . 
Such operators have dimension 15 or higher and have suppression factors of at least A~ n . 
These will lead to "triple nucleon decay" processes where three nucleons in a heavy nucleus 
undergo collective decays leading to processes such as pnn — > e + 7r°. We estimate the rates 
for such decay in Section 4 and find that A can be as low as 10 2 GeV. 
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Table 2: Charge assignment under B, L, Y and U(V)2Y-B+3L which contains the Zq of Tabled 



3 Embedding Z Q in l\ + Li + Lj — 2L k 

It is interesting to see if the Zq symmetry of Tabled can be realized as an unbroken subgroup 
of a gauged U(l) symmetry. Although the Zq is a subgroup of the U(1)(2y-b+3L), this U(l) 
would be anomalous without enlarging the particle content. We have found a simple and 
economic embedding of Z e into a U(l) gauge symmetry associated with 1^ + Lj + Lj — 2L&. 
Here Lj is the zth family lepton number and i ^ j ^ k. No new particles are needed to 
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cancel gauge anomalies. With the inclusion of righthanded neutrinos I R = Y — (B — L)/2 is 
an anomaly-free symmetry. Lj + Lj — 2L&, which corresponds to the A§ generator acting in 
the leptonic SU(3) family space, is also anomaly-free. 
The charges of the SM particles under this U(l) are 

Q t = (0,0,0), uf= (-1,-1,-1), d { c = (1,1,1), 

£, = (-4,2,2), ei c = (5,-1,-1), vf = (3,-3,-3) , H = l. 

This charge assignment allows all quark masses and mixings as well as charged lepton masses. 
When the U(l) symmetry breaks spontaneously down to Z 6 by the vacuum expectation value 
of a SM singlet scalar field <fi with a charge of 6, realistic neutrino masses and mixings are 
also induced. The relevant lagrangian for the righthanded neutrino Majorana masses is 

£„c = M 12 v{v c 2 + M 13 z> 3 c + u^ c 3 <f) + v c 2 v c 2 (j) + utul4>* + . (5) 

After integrating out the heavy righthanded neutrinos we obtain the following AL = 2 
effective operators: 

£al= 2 = j{iihHH + i&HH + t x t x HHe + £ 2 £ 2 HHe* + £ 2 £ 3 HHe* + £ 3 £ 3 HHe*). (6) 

Here A ~ M\ 2 ~ M\ 3 is the scale of L-violation and we have defined e = {(f)) /A. For e <C 1, 
this lagrangian leads to the inverted mass hierarchy pattern for the neutrinos which is well 
consistent with the current neutrino oscillation data. This neutrino mass mixing pattern is 
analogous to the one obtained from L e — L M — L T symmetry ^3]. However, here the U(l) is 
a true gauge symmetry. 

We have also investigated other possible U(l) origin of the Z 6 symmetry and found the 
I ft + Li + Lj — 2Lk combination to be essentially unique. To see this, let us assign a general 
U(l)x charge for zth generation of the SM fermions consistent with the Zq symmetry as 

{Qi, u% d% 4 e% vl} = {Qmf\ 5 + 6m?, 1 + 6m?, 2 + 6m?, 5 + 6m?, 3 + 6m?} 

where m? are all integers. The Higgs field has a charge H = 1 + 6m . If we impose the 
invariance of the Yukawa couplings of the charged fermions and Dirac neutrinos for each 



;eneration, the anomaly coefficients from ith generation become 
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(7) 

The coefficient for the mixed gravitational anomaly for each generation is zero. From Eq. 
©, it follows that A 2 = J2i A ?su<2) L fxu<i) x = + 9m ? + 3m 2 } ) = can be satisfied 



only when all three generation contributions are included. Once Ai = is satisfied, all 
other anomaly coefficients will automatically vanish. A 2 = can be rewritten in a familiar 
form as 3 J2i Q% + J2i ^ = 0- Thus we see that any U (1) symmetry satisfying this condition 
and consistent with the Z 6 charge assignment can be a possible source of Zg. If the Qi are 
different for different generations, quark mixings cannot be generated without additional 
particles. By making a shift proportional to hypercharge, we can set Qi = for all i. Two 
obvious solutions to J2i ^» = are ^ = (1) 1; — 2) and £; = (1, —1, 0). The latter one does not 
reproduce the Zq charge assignment while the former one does, which is our solution when 
1^ is added to it. 

A related B — 3L T has been discussed in Ref. [TH]. This is the same as B — L plus 
L e + L M — 2L T . In Ref , only one righthanded neutrino v c T is introduced so the seesaw 
mechanism applies only for one light neutrino. The other two neutrinos receive small masses 
from radiative correction. In our model, since there are three righthanded neutrinos, all the 
neutrino masses arise from the conventional seesaw mechanism. 



4 AB = 3 Operators and Triple Nucleon Decay 

We now list the lowest dimensional (d=15) AB = 3 effective operators which are consistent 
with the Zq symmetry. Imposing gauge invariance and Lorentz invariance, we find them to 
be: 



w c4 d c5 e c , w c2 J c V, Qu c3 d c5 £, Qu c2 d c6 I, Q 2 u c3 d c4 e c , 
Q 2 u c d c6 e c , Q 3 u c2 d c4 £, Q 3 u c d c5 £, Q 4 u c2 d c3 e c , Q A u c d c % c , 
Q A d cb e c , Q 5 u c d c3 £, Q 5 d c4 £, Q 6 u c d c2 e c , Q 7 d c2 £, Q 8 d c e c . (8) 

Here Lorentz, gauge and flavor indices are suppressed. These operators can lead to "triple 
nucleon decay". The dominant processes are 

PPP — > e + + 7T + + 7T + 

ppn —>■ e + + 7r + 
pnn — > e + + 7T° 

nnn — » v + 7r° . (9) 

Tritium ( 3 H) and Helium-3 ( 3 He) are examples of three-nucleon systems in nature. These 
nuclei are unstable and undergo /5-decay with relatively short lifetime. In the presence of 
operators of Eq. (jHJ), 3 H — > e + + tc° and 3 He — > e + + rr + decays can occur. However, there 
is no stringent experimental limits arising from these nuclei. So we focus on triple-nucleon 
decay in Oxygen nucleus where there are experimental constraints from water detectors. To 
estimate the decay lifetime we need to first convert the nine-quark operators of Eq. (jHJ) into 
three-nucleon operators and subsequently into Oxygen nucleus. 

We choose a specific operator Q 5 d c £/A u as an example to study the process pnn — > 
e + + 7r° triple nucleon decay process. This induces the effective three-nucleon operator in 
Oxygen nucleus 

Q 5 d c4 £ f3 3 {l + D + F) 



A 11 yftUK 



ii 



■(irnnpe) , (10) 
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where (3 ~ 0.014 GeV 3 is the matrix element to convert three quarks into a nucleon [To] . 
F ~ 0.47, D ~ 0.80 are chiral lagrangian factors, and = 139 MeV is the pion decay 
constant. 

We now estimate the wave-function overlap factor for three nucleons inside Oxygen nu- 
cleus to find each other. This is based on a crude free fermi gas model where the nucleons 
are treated as free particles inside an infinite potential well. A single nucleon wave function 
is given by ip m (x) = a/2/V sin(m7rx/r), where r is the size of nucleus and m is the energy 
level. Incorporating isospin and Pauli exclusion principle, the highest energy level which 
corresponds to m = 4 is found to have 2 protons and 2 neutrons. We assume the highest 
level has the most probability to form a Tritium-like "bound state" of three nucleons. The 
probability for three nucleons in Oxygen nucleus to overlap in a range the size of Tritium 
nucleus is 

P ~ A - f'^d m d d m (sin ( sin ( 4 ™t) sin ( *™) V ~ 0.0253 , 



3/n V r ' \ r J 



IT ' \ \ T* J \ f / \ f 



(11) 

where \ is the ratio between the radii of Tritium nucleus and Oxygen nucleus, since 

R oc A 1 / 3 (A is the atomic number). So the effective baryon number violating operator of 
Eq. (fTUj) becomes 

P8 3 

( 3 Hne) . (12) 



The triple nucleon decay lifetime can then be estimated to be 

167T^ 2 A 2 ^ 6 

r ~ pWm^ ' (13) 

By putting the current limit on proton lifetime of 3 x 10 33 yrs, we obtain: 

A ~ 10 2 GeV . (14) 

Thus we see the Zq symmetry ensures the stability of the nucleon. 

To test our crude model of nuclear transition, we have also evaluated the double nucleon 
decay rate within the same approach and found our results to be consistent with other more 
detailed evaluations [T7j . 

5 Realistic Neutrino Masses with Low Threshold 

If the threshold of new physics is low, but is still much larger than a TeV, small neutrino 
masses can be generated by the seesaw mechanism as discussed in Section 3 through the 
effective operators ££HH/A. If the threshold is as low as a few TeV, the induced neutrino 
mass will be too large. Here we show a mechanism by which such operators can be suppressed 
by making use of a discrete symmetry (with N odd) surviving to low scale. This Z^ has 
a natural embedding in B — L gauge symmetry. 
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Consider the following effective operators in the low energy lagrangian: 

CDUHH- + -. (15) 

Here S is a singlet field which has charge (1, 3) under Zn x Zq while £ has charge (—3, 2). 
(The Z§ charges of SM particles are as listed in Tabled) The first term in Eq. (fTKj) respects 
a (7(1) symmetry while the second term reduces this to Z 6 x Z^. If S develops a VEV 
of order 10 2 GeV, realistic neutrino mass can arise even when A is low. For example, if 
A = 10 TeV and S = 10 2 GeV, the neutrino mass is of order v 2 (S) G /A 7 ~ 0.4 eV, which is 
consistent with the mass scale suggested by the atmospheric neutrino oscillation data. 

Two explicit examples of the Zn symmetry with N = 5 and 7 are shown in Table El 
These Zn symmetries are free from gauge anomalies. In the Z5 example, the crossed anomaly 
coefficients for SU(3)c and SU (2) l are 5N g and 5N g /2 respectively showing that Z 5 is indeed 
anomaly-free. For Z?, these coefficients are 7N g and 7N g /2, so it is also anomaly-free. 
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Table 3: Zjv charge assignment for N = 5 and 7. 



It is interesting to ask if the Z^ can be embedded into a gauged U(l) symmetry. A 
simple possibility we have found is to embed this Z^ into the anomalous U(1)a symmetry of 
string origin with the anomalies cancelled by the Green-Schwarz mechanism [TH]. Consider 
U(1)b-l without the right handed neutrinos but with the inclusion of vector-like fermions 
which have the quantum numbers of 5(3) and 5(2) under SU(5) x U{1)a- This U{1)a is 
anomaly-free by virtue of the Green-Schwarz mechanism. When this U(1)a breaks down to 
Z5, the extra particles get heavy mass and are removed from the low energy theory which is 
the Zq x Z 5 model. 

Without the second term in Eq. (fTKjl . the phase of S field will be massless upon spon- 
taneous symmetry breaking. This Majoron field ^H] would however acquire a mass from 
the second term of Eq. ()15|1 . In the Z$ x Z^ model, the mass of the Majoron is of order 
(S) 7 /A 6 ~ 100 keV. In the Z 6 x Z 7 model, the Majoron mass is of order (S) u /A 10 ~ 10 eV. 
Such a Majoron with a mass of either 100 keV or 10 eV is fully consistent with constraints 
from early universe cosmology [20 . The interaction term ££HHS 6 /A 7 induces the Majoron 
decay S — > vv with a Yukawa coupling Ys-* w = 6m^/ (S) ~ 10~ n . The decay rate of the 
Majoron can be estimated to be 

r = Y ^ mS ~ 10- 23 m 5 • (16) 

This corresponds to a Majoron lifetime of r ~ 10 sec for the Z§ x Z5 model and r ~ 
10 5 sec for the Z§ x Z7 model. Such a Majoron can modify the big-bang nucleosynthesis 
processes. However the modification is not significant since the Majoron will decouple before 
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the electro-weak phase transition. Its contribution to the expansion rate is equivalent to 
that of 0.047 x 4/7 ~ 0.027 light neutrino species [21J. This extra contribution is well within 
observational uncertainties. 



6 Conclusions 

We have shown in this paper that the SM lagrangian including small neutrino masses and 
with three generations has an anomaly-free Z§ discrete symmetry. This Z§ can act as gauged 
baryon parity and stabilize the nucleon even when the new physics threshold is below the 
TeV scale. The Z 6 has a natural embedding in l\ + Li + Lj — 2Lk gauge symmetry. We have 
also shown that realistic neutrino masses can be obtained through an additional anomaly-free 
Z 5 or Z 7 symmetry which may have its origin in the anomalous U(1)a symmetry. 

Our framework can be readily extended to incorporate supersymmetry. The Z 6 charge 
assignment is still as in Table H The down-type Higgs H has a Z 6 charge of 5 and the 
Grassmann variable 6 has charge 3 or ^U]. The Higgsino contribution to the [SU(2) j,fxZ 6 
anomaly is 3 and the gaugino contribution to [SU(2) L } 2 x Z 6 and [SU(3)c] 2 x Z 6 anomalies 
are 6 and 9 when 9 charge is not zero. Since these contributions are all multiples of 6/2, 
the Zq symmetry remains anomaly-free in supersymmetric case. This Z§ forbids all AB = 1 
and AB = 2 nucleon decay as in the non-supersymmetric case. 
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